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Highly confined vectorial electromagnetic 
field distributions represent an excellent 
tool for detailed studies in nano-optics and 
high resolution microscopy, such as nonlinear 
microscopy^, advanced fluorescence imaging^^ 
or nanoplasmonics^'^. Such field distributions 
can be generated, for instance, by tight focussing 
of polarized light beams^^^. To guarantee high 
quality and resolution in the investigation of 
objects with sub-wavelength dimensions, the pre- 
cise knowledge of the spatial distribution of the 
exciting vectorial field is of utmost importance. 
Full-field reconstruction methods presented so 
far involved, for instance, complex near-field 
techniques^^'^^. Here, we demonstrate a simple 
and straight-forward to implement measurement 
scheme and reconstruction algorithm based 
on the scattering signal of a single spherical 
nanoparticle as a field-probe. We are able to 
reconstruct the amplitudes of the individual focal 
field components as well as their relative phase 
distributions with sub- wave length resolution 
from a single scan measurement without the 
need for polarization analysis of the scattered 
light. This scheme can help to improve modern 
microscopy and nanoscopy techniques. 

In the optical analysis of sub- wavelength objects such 
as cellular structures'^, plasmonic particles^'^'^^ or sin- 
gle spins^, nano-optical tools including highly resolving 
microscopy techniques are used. Because such methods 
utilise complex and highly confined vector fields, the ex- 
act knowledge of the corresponding spatial field distribu- 
tions is crucial. In the last decades, several techniques 
have been proposed to map these focal fields, such as 
using metal knife edges'^'^^ to probe the total electric 
energy density distribution, fiuorescence molecules'^, ta- 
pered fibres'^' '^ or tip-based methods'^ to image specific 
field orientations, or near-field scanning optical micro- 
scope (NSOM) techniques to extract amplitude and even 
phase information''. These NSOM-based methods re- 
quire complex measurement and detection schemes and 
calibration procedures to allow for amplitude and phase 



mapping of individual field components. As an alter- 
native approach for measuring phase information also a 
single particle scattering scheme was proposed recently^^ , 
where the authors show that Mie-scattering can distin- 
guish the topological charge of vortex-beams. 

We now demonstrate a precise and easily imple- 
mentable field reconstruction technique for highly con- 
fined field distributions created by arbitrary focusing sys- 
tems, based on, what we call, Mie-scattering nanoint- 
erferometry The basic concept of this reconstruction 
method for highly confined focal field distributions can 
be understood as follows. We use a metallic nanosphere 
on a glass substrate as local field sensor and scan it step- 
wise through the focal field distribution under investiga- 
tion. If the particle size is chosen small enough, it will 
be excited by the local electric field E(ro) only. Its re- 
sponse can be described, in a first approximation, by an 
electric dipole. Hence, for a particle on a substrate differ- 
ently oriented electric dipole moments will be excited in 
the particle for different positions in the field distribution 
under investigation depending on the corresponding lo- 
cal field orientation, and resulting in different scattering 
patterns^^ (see Fig. 1(a)). Here we show, that the initial 
local field and therefore the three-dimensional focal field 
distribution can be determined accurately in amplitudes 
and relative phases of the individual field components by 
collecting the transmitted light angularly resolved^' and 
exploiting the interference signal between incoming and 
scattered field. The interference and, hence, the phase 
information is preserved herein by effectively changing 
the observation direction. 

To analytically describe the full scattering process, 
the unknown focal field is expanded into electromagnetic 
multipoles. With this choice of a basis system, only the 
lower expansion orders have to be taken into account for 
highly confined fields (see supplementary section S.6 and 
23 and 24). The focal electric field distribution can thus 
be expressed as 
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where N^^ and M^^ are regular vector spherical har- 
monics representing differently oriented electric and mag- 
netic multipoles^^'^^ expanded around the geometrical fo- 
cus. This transfers the full information of the unknown 




the known geometry and refractive indices as: 



FIG. 1: Schematics of the scattering process, a Far-field 

emission pattern of differently oriented electric dipoles 

on an interface with refractive index n = 1.5. The 

dipole intensity emitted into a limited solid angle 

(dashed lines) shows a strong dependence on the dipole 

direction, b The interface additionally leads to a 

coupling of the reflected field from the substrate 
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Ein/s and the scattered field off the spherical 
scatterer Eo. 



focal field from the complex-valued electric field com- 
ponents Ein,i(r) at each point in the focal plane to the 
complex- valued multipole expansion coefficients Amn and 
Bmn- This representation allows to relate distinct far- 
field patterns with the focal field distribution under inves- 
tigation, where electric dipoles can be associated with the 
local electric field Ein(ro) at the expansion point tq, elec- 
tric quadrupoles and magnetic dipoles with its gradient, 
and higher order multipoles with higher order moments 
of the field. In addition, it permits a simple treatment of 
scattering problems via the T-matrix approach^^. 

The focal electric field distribution can thus either be 
reconstructed by determining all multipole coefficients of 
the focal field at one single point or by scanning the dipole 
contributions in the focal plane and relating them to one 
common expansion point via the translation theorem for 
vector spherical harmonics^^. To achieve the in both 
cases necessary high position accuracy and repeatability, 
the probe is immobilized on a substrate in our exper- 
iments, which introduces an interface in the scattering 
problem (see Fig. 1). The resulting transmitted field Et 
can then be expressed as the incoming field Ein,t and a 
scattering field term Es,t, incorporating the full geometry 
of the system: 



Et(r) = Ein,t(r)+Es,t(r). 



(2) 



For a simplification of the description of the interface, we 
transform the basis functions in the forward hemisphere 
to N^n,t and M^n,t, leaving the expansion coefficients 
Amn and B^^ unaltered and thus keeping the low or- 
der multipole expansion (see supplementary section S.l). 
The scattered electric field Eg^t is then related to the 
incoming field Ein,t by describing the interaction of the 
nanoprobe and the substrate with the incoming light field 
by an effective scattering matrix Teff , depending only on 



E,,= (l-TLlf))-^T(l + LW) 



Ein,t — TeffEin,t, 



(3) 



with T being the scattering matrix of the probe parti- 
cle (including information about the particle shape, size 
and its optical properties) and L)^' ^ being reflection op- 
erators of the substrate (see Fig. 1(b) and supplementary 
section S.l). 
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FIG. 2: Sketch of the experimental implementation of 

the reconstruction scheme, a A spherical metallic 

nanoparticle adhered on a glass substrate is scanned 

through the focal field distribution under investigation. 

The scattered and transmitted intensity is consequently 

collected with a variable solid angle (Oa^^Pa)- b 

Experimental setup: an arbitrarily structured input 

beam is focused by a high numerical aperture 

microscope objective (NA = 0.9). The nanoprobe can 

be precisely scanned through the focal plane by means 

of a 3D-piezo stage. The transmitted light is then 

collected by an immersion-type microscope objective 

(NA = 1.3) and its back- focal-plane is imaged onto a 

C CD-chip. The reflected light can be measured in 

addition to ensure the exact focus position. 



To relate the unknown input field to an experimentally 
measurable quantity, the resulting power transmitted at 
an angle {9^(f>) (see Fig. 2(a)) can then be expressed by 

Pt(^, V) = Pi„(0, V) + Ps(^, ^) + Pext(^, ^), (4) 

following the classical Mie scattering problem^^with 
Pin = ^Re [Er„,t X Hi„,t] , Ps = ^Re [E* ^ x H,,t] , 
Pext = ^Re [Er„,t X H,,t + E* , X Hi„,t] 

as the incoming, scattered and extinct power measured 
in transmission. Here, the magnetic field components are 
determined from the electric field components in the far- 
field using the plane wave spectrum of the transmitted 
multipoles (see supplementary section S.l). The interfer- 
ence term Pext, which depends on both and (/9, not only 
allows for the extraction of amplitude information of the 
multipole expansion coefficients, but also the phase rela- 
tion between them. It can be shown, that for an unam- 
biguous reconstruction of both, amplitudes and relative 
phases of the individual electric field components in the 
focal plane the combination of the following collection 
schemes is sufficient. First, the transmitted power Pt 
has to be recorded and integrated in a given solid angle 
around the optical axis for each position of the particle 
in the focal plane (see Fig. 3(a) top and center). For 
that purpose, an integration in from to Oa and in 
Lp from to 27r is performed. By this integration over 
the whole range in (/p, some phase information is lost. 
Thus, to regain the full information about the relative 
phases, an effective break of the cylindrical symmetry 
of the collection system (objective) is introduced. This 
can be realized by secondly choosing a collection sector, 
hence performing a 6>-integration from to Oa and a ip 
integration from ipi to ip2 (with |(/:'i — ^2! < 27r; see Fig. 
3(a) bottom). With this choice, the light emitted into a 
certain direction relative to the optical axis is analysed. 
Both schemes can be realised easily in the experiment 
by imaging the transmitted power angularly resolved in 
the back focal plane of a high numerical aperture (NA) 
collection lens for each position of the nanoprobe rela- 
tive to the beam. The integration ranges can then be 
chosen freely, with bound by the maximum collection 
angle ^max of the employed collection lens. Additionally 
collected angular ranges can furthermore reduce the in- 
ffuence of noise in the system (see supplementary section 
S.2). With the known scattering matrix Teff, the powers 
transmitted into different angular ranges for each posi- 
tion of the nanoprobe relative to the input beam thus 
establish a system of positive quadratic forms of the un- 
known multipole expansion coefficients. The inversion 
of this equation system then leads to the vectorial focal 
electric field distribution. 

For the experimental demonstration of the introduced 
nanointerferometric vectorial field reconstruction scheme 
a custom-made scanning setup (similar to Ref. 11) is 
used, which represents a scattering system that is easily 



applicable to different research areas. The highly con- 
fined focal field distribution under investigation is cre- 
ated by sending a collimated input beam with a wave- 
length of 530 nm into a high NA objective. The probe 
system which we use for scanning consists of a spherical 
gold nanoparticle with a diameter of 82 nm on a glass 
substrate (see Fig. 2). The diameter of the particle was 
chosen such, that it is small enough to exhibit predomi- 
nately a dipolar response but still large enough to guar- 
antee a large scattering cross-section for a high signal to 
noise ratio. Furthermore, the spherical symmetry of the 
particle simplifies the scattering matrix T. The parti- 
cle is scanned through the focal field distribution and at 
each point the light transmitted through the substrate 
and scattered off the nanoprobe in the forward direction 
is collected by an immersion-type microscope objective 
(see Fig. 2(b)). The back focal plane of the collection 
objective is then imaged onto a C CD-camera, provid- 
ing access to the k-spectrum of the transmitted light. 
Integrating this angularly resolved transmitted intensity 
over different angular ranges {Oa^^a) (see Fig. 2(a)), en- 
ables the reconstruction of the vectorial focal field dis- 
tribution as discussed above. The collection objective is 
hereby assumed to only introduce negligible additional 
errors to the collected intensity distribution, which can 
be achieved with modern high quality oil-immersion ob- 
jectives. It is worth noting here, that the aforementioned 
scheme is also applicable to refiection-type setups and 
probes embedded in a homogeneous medium. 

By choosing a radially polarized doughnut beam as 
exemplary input field, a three-dimensional vectorial field 
distribution is generated under tight focusing conditions, 
exhibiting an on- axis longitudinal field component^' ^^ 
and off-axis transverse field components (see Fig. 3(c), 
calculated via vectorial diffraction theory^). Such a 
tightly focused cylindrical vector beam finds applications 
in several fields of nano-optics and imaging^^'^^. 

In Fig. 3(a), the experimental results are shown. For 
every position of the nanoparticle (scan step-size: 25 nm; 
see supplementary section S.5) a single image of the back 
focal plane is recorded (see Fig. 3(a) left column for one 
position of the particle and supplementary section S.3). 
From this camera data, two-dimensional scan images are 
derived by plotting the transmitted intensity integrated 
over the corresponding angular range in the back focal 
plane image for each probe position (see Fig. 3(a) right 
column). Here, we choose three different angular ranges 
for integration in the measured back focal plane images 
to reduce the inffuence of experimental noise (see supple- 
mentary section S.2). We choose two solid angles which 
correspond to two full numerical apertures of 0.9 and 
0.4 and a sector of (/:?o = 1 rad at an NA of 0.9. The 
non-rotational symmetric collection angle preserves the 
interference information (see Fig. 3(a) bottom). This 
signal is highly sensitive to the actual optical proper- 
ties of the nanoprobe, which were determined experi- 
mentally for the selected nanoprobe and wavelength to 
cauNP = —3.0 + 2.12 (see supplementary section S.4). 
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FIG. 3: Experimental results and theoretical comparison for a radially polarized vector beam, a Image of the back 
focal plane of the collection objective with an NA of 1.3 for one position of the nanoprobe relative to the focal field. 

The measured intensity scan-images in transmission for a wavelength of 530 nm correspond to three different 
collecting solid angles with a full aperture of NA 0.9 and 0.4 as well as an azimuth angle oi ipa — ^ rad for an NA of 

0.9, all derived from the same measured back focal plane images for different probe positions. The intensity is 

normalized to the total intensity of the input beam, b Squared electric field components |£^i| and relative phases <l>i 

in the focal plane reconstructed from the measured intensity distributions up to a multipole order of n = 8. c Energy 

density distribution for the three field components in the focal plane of the same beam calculated via vectorial 

diffraction theory. The insets show the calculated phase distribution in the same scale and colormap as in (b). 



The electric energy densities and phases of the focal field 
components reconstructed from these intensity distribu- 
tions (see Fig. 3(b)) show a very good overlap with the 
calculated field components shown in Fig. 3(c). Small 
deviations from theory are visible in the profiles recon- 
structed form the experimental data. We attribute these 
differences to small aberrations of the focussing system. 
Thus, deviations from the ideally expected focal field dis- 
tribution can be reconstructed with sub-wavelength pre- 



In summary, we have shown an easily applicable re- 
construction scheme to determine the full vectorial am- 
plitude and relative phase distributions of highly con- 
fined electromagnetic fields. The technique relies on the 
nanointerferometry between the input field and the field 
scattered off a nanoprobe as well as an angularly re- 
solved measurement of the resulting far-field intensity. 
By adapting the described scheme, also a reconstruction 
of field distributions for other focusing systems can be 
realised, including the near-field distribution of NSOM- 
tips. 



METHODS 



Experimental Setup 

A tunable light source (NKTphotonics SuperK Power- 
Plus with an AA Opto-Electronic MDSnC-TN acousto- 
optical tunable filter) is providing a linearly polarized 
Gaussian beam at a wavelength of 530 nm, which is op- 
tionally converted into a radially or azimuthally polar- 
ized doughnut beam by a liquid crystal radial polariza- 
tion converter (ARCoptix). The beam is then guided via 
four mirrors top-down onto a microscope objective with a 
NA of 0.9 (Leica HCX PL FL lOOx/0.90 POL 0/D). The 
beam radius (1/e^) of the doughnut beam is 1.64mm, 
while the diameter of the entrance aperture of the mi- 
croscope objective is 3.6 mm. A single spherical gold 
nanoparticle with a diameter of 82 nm sitting on a glass 
cover slip is scanned through the focal plane by a high 
precision 3D piezo table (PI P-527) and the transmitted 
light is collected via an oil immersion objective with a 
NA of 1.3 (Leica HCX PL FLUOTAR lOOx/1.30 OIL). 
In addition, the reflected light is collected by the focusing 
objective. The angular distributions of the transmitted 
and optionally the reflected light are detected via imag- 
ing the back focal plane of the corresponding objectives 



on CCD-cameras (The Imaging Source DMK 23G618). 



Sample preparation 

The spherical gold nanoparticle was fabricated at the 
Max Planck Institute for Polymer Research, Mainz, Ger- 
many by laser induced melting of a commercial col- 
loidal gold solution. The resulting solution of spherical 
nanoparticles was drop-coated on a microscope cover slip 
with a thickness of 170 pm and prestructured gold mark- 
ers to reproducibly measure one single nanoparticle. 
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SUPPLEMENTARY MATERIALS 



basis vectors assume the form 



This supplementary material describes in more de- 
tail the developed reconstruction algorithm. In addi- 
tion, we show the characterisation results of the employed 
nanoprobe as well as the reconstruction of a tightly fo- 
cused azimuthally polarised doughnut beam. Further- 
more, the theoretical and experimental sources of errors 
are discussed. 
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S.l. Detailed description of the influence of the 
substrate 



Here, k^^ is the wave vector of a reflected plane wave, 
e^^^ are unitary vectors of the reflected coordinate system 
and Rsp is the radius of the nanoprobe. To account for 
the interaction of the reflected fleld with the nanoprobe 
above the interface, the reflected multipoles (6) have to 
be reexpanded into the non-reflected ones. Thus we write 



Describing the scattering process off the nanoprobe in 
the basis of vector spherical harmonics (VSHs), we follow 
the notation of Tsang et al.^. To simplify the theoretical 
description of the energy transfer through the interface 
and the interaction of the scattered flelds with the latter, 
we transfer the effect of the interface to the basis func- 
tions M.lin and N^^n, leaving the expansion coefficients 
of the initial focal fleld unaltered. For this purposes we 
need a plane wave representation of the flrst and third 
type VSHs, which results in^ 
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where jmn are multipole dependent prefactors, k is the 
wave vector of the respective plane wave and X^nn and 

("2) 

X^n correspond to VSHs defined analogue to Ref. 1. 
The index j distinguishes the incoming and scattered 
fields, using VSHs of first and third type. The integra- 
tion region of these two types differs due to the inclusion 
of evanescent field terms in the scattered fields and is 
ni = (f e (0, 27r) , 6> G (0, tt), and Qs = ^ ^ (0, 27r) , 6> G 
(0, 7r/2 — ioo) if an observation point is in the upper hemi- 
sphere OT Qs = (f e (0, 27r) , e (7r/2 + ioo, tt) if it is in 
the lower hemisphere. 

The effect of the interface can now be introduced by 
applying Fresnel equations to each plane wave component 
in Eq. (5). We here use the Fresnel transmission t^ (6) 
and reflection r^ {0) coefficients from Ref. 2 (where rj = 1 
for or p-polarised components and r] = 2 for (f or s- 
polarised components) . As a result of the reffection the 
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The elements of the reffection operators L)^' can then 
be readily obtained from Eq. (7). 

Concerning the initial and scattered fields transmitted 
through the interface, the free space multipoles can be 



similarly expressed as 
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Thus, when an electromagnetic field represented as a vec- 

Brr 



tor 



A.n 



containing expansion coefficients of the free- 



space multipoles is transmitted through a planar inter- 
face, the electromagnetic multipoles are replaced by their 
transmitted (10) counterparts but the expansion coeffi- 
cients of the transmitted electric field remain unchanged. 
This fact enables us to write the density of the Poynting 
vector in the hemisphere of the substrate in a matrix 
form as 
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(11) 



(T) (T) (T) 

where the matrices w^^ , Wg and We contain scalar 
products of the plane wave spectrum functions of the 
transmitted multipoles. The energies transmitted to a 
solid angle (Oa^^fa) can then be obtained by performing 
an integration of those matrices over the desired solid 
angle and applying (11). 



S.2. Mathematical background of the 
reconstruction algorithm 

From a mathematical point of view the density of the 
Poynting vector as represented in Eqs. (11) is a posi- 
tively defined quadratic form, which has a unique solu- 
tion. However, when an integration is performed over a 
given solid angle, the rank of the matrix may decrease 
(for example, if cp = (0,27r)). Basically this means, that 
for some combinations of m and n the crossproducts of 
the multipole amplitudes (like AijA*,-, or AijB^,j,) do 
not infiuence the experimentally observed transmission. 
Therefore, in order to preserve the inter ferometric infor- 
mation, the integration region has to be chosen carefully, 
so that the rank of the matrix does not decrease. Then, 
we end up with a number of quadratic equations, which 
can be solved using a variety of different techniques. 

The most straightforward inversion approach is based 
on a relinearisation of quadratic equations^, where one 
introduces new variables Fi = AiAj (here Ai is either 
the real or imaginary part of the multipole amplitude). 



which are linearly related to the transmission at each 
scan position. In an ideal case, the unknowns Fi are 
trivially retrieved and thus the unknown multipole am- 
plitudes are obtained. However, in an experimental sit- 
uation, one always has to account for some systematic 
and statistical perturbations, therefore a straightforward 
matrix inversion of the form Ax — y = ^ may fail. De- 
pending on the strength of the stochastic term ^ relative 
to the real signal, some quadratic forms may become sin- 
gular in the sense, that some crossproducts for a given 
integration region inffuence the transmission signal less 
than it is inffuenced by the stochastic term. This can be 
counteracted by a careful choice of additional integration 
regions, where that crossproduct contributes stronger to 
the transmission. 



S.3. Sampling effect in the experimental 
CCD-images 

To measure the energies transmitted to certain solid 
angles experimentally, the back focal plane of the micro- 
scope objective collecting the transmitted and forward 
scattered light is imaged onto a C CD-camera. This in- 
troduces a sampling effect on the recorded back focal 
plane images, which has to be corrected when comparing 
the resulting experimental and theoretical transmitted 
energies. The sampling results from the ffnite pixel size 
of the CCD-camera being related to ffnite integration an- 
gles in the spherical coordinate system of the transmitted 
multipole components. Thus, the effective integration 
area of each pixel gets transformed from dF = Ax Ay 
to dF' = cos6A6A(p (see Ref. 4). Measuring the elec- 
tric energy density lEl^^^^j^, the energy transmitted at 
a certain collection angle is then |E| = cos^ 6> lEl^^^^j^. 
The sampling effect was additionally checked for a ra- 
dially polarised input beam of known size by rigorously 
calculating the resulting ffelds through the focusing and 
collection objective and comparing the theoretically ex- 
pected energy density with a CCD-image of the trans- 
mitted beam (see Fig. S 1). 



S.4. Characterisation of the used nanoprobe 

The precise knowledge of the shape and optical prop- 
erties of the used nanoprobe govern the accuracy of the 
proposed reconstruction algorithm. Since the employed 
gold nanoparticle is fabricated via laser-induced melt- 
ing, we result in a spherical nanoprobe of 82 nm diam- 
eter, measured via SEM. To account for possible changes 
in the material conformation, we measured the effective 
refractive index of the nanoprobe via Mie-scattering^ . 
Thus, the nanoprobe is embedded in a homogeneous 
medium via index-matching immersion oil (n = 1.52) 
and the focussing objective is exchanged with a second 
oil-immersion objective with a numerical aperture (NA) 
of 1.3. By choosing a linearly polarised Gaussian input 




Figure SI: Influence of sampling: comparison between 
the theoretically expected beam cross-section in the 
back focal plane of the collection objective and the 
experimentally retrieved cross-section for a radially 

polarised input beam with wq/ f = 0.69, corrected by 
the sampling- factor cos^ 6. 
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Figure S2: Determination of e: Measured and fitted 

reflectance (top) and transmittance (bottom) spectrum 

for the used nanoprobe embedded into a homogeneous 

medium with n^ = 1.52. 



beam with a small beam width of wo/f = 0.35, the effec- 
tive NA of the system is reduced to 0.5. Positioning the 
particle in the focal spot of the system, an approximately 
homogeneous plane wave excitation is achieved over the 
whole particle. 

The normalised intensity collected in reflection and 
transmission is then linked to the scattering and absorp- 
tion cross-section of the particle via 
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l-Ao-(Q 
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':Qs 



(12) 



where Aq depends on the polarisability of the particle. 
The factors | and | follow from the maximum collection 
angle of ^max = 60°, where approximately one third of 
the scattered intensity of a transverse dipole is collected 
in forward direction as well as in backward direction. 

To retrieve the relative permittivity of the particle, the 
transmittance and reflectance of the embedded particle 
are measured for a spectral range of A = 475 — 710 nm. 
This spectral data is then used to fit a physical model of 
the particle's relative permittivity to the classical Mie- 
scattering solution of Qsca and Qabs- The utilised model 
is based on an extended Drude model with critical point 
description^ for the interband transitions and was shown 
to give excellent results for fitting the permittivity of pla- 
nar gold films. The fitted spectral response shows a very 
good overlap with the measured resonance spectra (see 
Fig. S 2) and results in 

e(A = 530 nm) = (-3.0 ± 0.2) + (2.1 ± 0.2)z 



S.5. Effect of scanning and influence of the step size 

When scanning the nanoprobe through the focal field, 
one has to consider the general effect of translating the 
probe as well as the influence of the chose step size. 
Thus, let us consider one scanning step by a new co- 
ordinate frame {R' ^0' ^cp')^ related to the old one by 
R = R' + Rq. Here, Rq = (i^o^^o^^o) is an arbitrary 
displacement, where the polar and azimuthal coordinate 
axes remain unchanged. The functions M^^, N^^ of 
the old coordinate frame can then be expressed as a 
sum of the functions MJ^^, N|J^^ of the new coordinate 
frame. The corresponding electric field is then expressed 
as Ei„ = Er=i E;:=-. A'^,N'^,+B'^,M'^,, where the ex- 
pansion coefficients A^^^^ and B^^^ can be calculated from 
the addition theorem of VSHs^. In a matrix representa- 
tion, this corresponds to 
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A 



(13) 



where Trij(Ro) are sub-matrices of the translation oper- 
ator Tr(Ro). 

To determine an optimized scan step Ap < A for the 
experiment, the elements of a translation matrix are usu- 
ally represented as the sums 
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(14) 



for the wavelength used in the reconstruction scheme. 
The resulting wavelength dependent permittivity addi- 
tionally reproduces the resonance wavelength of colloidal 
gold particle solutions with different diameters very well 
(tested with particles from nanoComposix, Inc.). 



where (fo denotes the azimuthal coordinate (direction) of 
the scan step, \u - n\ < p < u ^ n, a^^(jp) and h^^(p) 
are transition matrix elements^, ji{x) denotes the spher- 
ical Bessel functions of first kind, and P/ {x) stands for 
the associated Legendre polynomials. For a small steps 



Ap ^ 0, the spherical Bessel functions can be approxi- 
mated to jo ^ 1 + 0(Ap2), ji ^ A:Ap/3 + 0(A/)^), and 
for n > 2 to j^ ~ O(Ap^). Translating the nanoprobe 
by a small distance, the sums in Eq. (14) thus contain 
only two terms: p = with m = fi and p = 1 with 
\m — ii\ = 0,1. This implies that only electric dipoles 
and quadrupoles as well as magnetic dipoles have to 
be considered. Since the employed nanoprobe responds 
mainly to electric dipoles, the size of the step has to be 
chosen, so that a translated magnetic dipole or electric 
quadrupole results in a measurable change of the electric 
energy density in the focal plane. This statement can be 
mathematically formulated as Tr = 1 + /cApOi, where 1 
is a diagonal matrix and /cApOi represents small linear 
changes. The size of the step now has to be chosen, so 
that the change in the scattering energy /SWgca due to 
the change of the coordinate system is larger than the 
background noise VKnoise introduced in the experimental 
measurements. 

Our estimation shows, that for our nanoprobe system, 
a scan step of Ap = 25 nm is large enough to tolerate 
fluctuations in the energy of up to 8%, whereas scan steps 
of Ap = 10 nm are sufficient to tolerate fluctuations of 
up to 3%. 



S.6. Additional sources of uncertainties 

Due to the desired sub- wavelength resolvability of the 
focal field reconstruction, several experimental and nu- 
merical uncertainties have to be considered. From an 
experimental point of view, the exact knowledge of the 
optical axis of the collection objective as well as its op- 
tical properties concerning aberrations and polarisation 
dependent transmission coefficients have to be exactly 
known. A mispositioned optical axis corresponds to a 
rotation of the chosen solid angles and thus introduce 
a changed weighting of the field components. Pinpoint- 
ing the axis position in the back focal plane to better 
than ko/100 results thus in an accurate reconstruction of 
the electric field components up to a relative strength of 
10~^. Any strong aberration of the collection objective 
on the other hand can not be easily distinguished from 
aberrations in the reconstructed field and thus lead to 
ambiguities. We found that for current optimized high- 
NA oil immersion objectives this error source is negligi- 
ble. 

Concerning the nanoprobe, its exact optical material 
parameters and its geometrical shape have to be deter- 
mined accurately, as mentioned in S.4. A spherical probe 
is of distinct advantage here, making the size estimation 
via SEM and the material parameter determination via 
Mie-scattering routinely feasible. 

From a theoretical point of view, the maximum mul- 
tipole order of the focal field sets a limit to the lateral 
dimension of the focal field as well as the accuracy to re- 
construct small structures of the field. Due to the trans- 
lation of the probe through the focal field, the expansion 
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Figure S3: Relative power of a radially polarised focal 

field (with the same parameters as in the main 

manuscript) contained in the ffist a n^ax = 5, b 

^max = 8 and c nmax = 10 multipole orders when being 

translated through the focal plane. The lateral position 

corresponds to a relative shift of the beam axis with 

respect to the expansion axis of the multipoles. 



point of the scattered and extinct fields is shifted with 
respect to the optical axis. Thus, power is distributed to 
higher order multipoles, leading to a distorted scattering 
response when considering only n multipole orders. This 
effect can be seen in Fig. S 3, where the power contained 
in 5, 8 and 10 multipole orders is shown for a radially 
polarized focal field with the same parameters as in the 
main manuscript when the beam is shifted with respect to 
the central position. While 5 multipole orders still show a 
significant drop in power at the outer borders of the scan 
area, 8 multipoles give already a fairly good account of 
the beam power. With 10 multipole orders, the influence 
of the lateral shift is then negligible for the dimensions of 
the actual focal field. Considering tightly focused fields 
with focal spot diameters smaller than the wavelength, 
a multipole order of nmax = 8 thus shows accurate re- 
construction results within the focal spot (see also Refs. 
8 and 9). Additionally, the choice of the collected solid 
angles is crucial to achieve unambiguous results as shown 
in S.2. 



S.7. Reconstruction of further vectorial focal field 
distributions 

To conffim the capabilities of the introduced vecto- 
rial reconstruction scheme, we experimentally measured 
different tightly focused field configurations. Figure S 
4 shows the experimental results of an azimuthally po- 
larised input beam for three different ranges of integra- 
tion in the back focal plane of the collection objective 
(Fig. S 4(a)) as well as the reconstructed (Fig. S 4(b)) 
and theoretically calculated (Fig. S 4(c)) focal field dis- 
tributions. The reconstructed and numerically evaluated 
electric field components show an excellent overlap, with 
only a weak residual longitudinal electric field compo- 
nent stemming from either uncertainties in the numeri- 
cal reconstruction or noise of the measurement. From the 
integrated intensity collected in a solid angle of 1 rad, de- 
tails of the interferometric information down to a lateral 
variation of 100 nm or A/5 can be resolved. 
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Figure S4: Results for an azimuthally polarized vector beam, a Image of the back focal plane of the collection 

objective with an NA of 1.3 for one position of the nanoprobe relative to the focal field. The measured intensity 

scan-images in transmission for a wavelength of 530 nm correspond to three different integration solid angles with a 

full aperture for an NA of 0.9 and 0.4 as well as an azimuth angle of cpa = I rad for an NA of 0.9. The intensity is 

normalized to the total intensity of the input beam, b Squared electric field components |£^i| and relative phase $i 

in the focal plane reconstructed from the measured intensity distributions, c Energy density distribution for the 

electric field components in the focal plane of the same beam calculated via vectorial diffraction integrals. The insets 

show the calculated phase distribution in the same scale and colormap as in (b). 
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